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An approximate method is proposed to solve position dependent mass Schrodinger equation. The 

IT) ' 

procedure suggested here leads to the solution of the PDM Schrodinger equation without trans- 
forming the potential function to the mass space or vice verse. The method based on asymptotic 
Taylor expansion of the function, produces an approximate analytical expression for eigenfunction 

-i— > ■ 

and numerical results for eigenvalues of the PDM Schrodinger equation. The results show that 



PDM and constant mass Schrodinger equations are not isospectral. The calculations are carried 
out with the aid of a computer system of symbolic or numerical calculation by constructing a 



simple algorithm. 
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INTRODUCTION 



Quantum mechanical systems with a position dependent mass (PDM) generate interest 
for its relevance and importance in describing the physics of many microstructures of current 
interest, understanding transport phenomena in compositionally graded crystals, designing 
modern fabrication of nano devices such as quantum dots, wires and wells, developing the- 
oretical models for effective interactions in nuclear physics, neutron stars, liquid crystals, 
metal clusters [l- 15]. These applications have stimulated a naturally renewed interest in 



the solution of PDM quantum mechanica 
Schrodinger equation, Dirac equation 16 



Hamiltonians. Recently, solution of the PDM 



191 ] and Klein-Gordon equation 2CH22j have re- 



ceived much attention. A number of authors have studied PDM Schrodinger equation within 
the framework of point canonical transformations 
super symmetric quantum-mechanical 



35 



23M26I]. Lie algebraic techniques [2 7H34|. 



431 ] . or other related techniques [44463]. 



In most applications of such methods, PDM Schrodinger equation has been transformed 
in the form of the constant mass Schrodinger equation by changing coordinate and wave func- 
tion. Obviously, this transformation generates isospectral potentials and exact solvability 
requirements result in constraints on the potential functions for the given mass distributions. 
In other words, a suitable transformation of coordinate and wave function becomes a bridge 
between constant mass and position dependent mass Schrodinger equation. As an example 
in a constant mass Schrodinger equation the choice of coordinate u = J® \Jm(x)dx and 
wave function ip(u) = [2m(x)] 1//4 tp(x) provides its transformation in the form of the PDM 
Schrodinger equation. In this case the potential is mass dependent; i. e. harmonic oscillator 
potential can be expressed as V = \muj 2 u 2 = \muj 2 ^ J* ^Jrn\x)dx^ and both constant and 
PDM Schrodinger equations have the same eigenvalues. The origin of such an isospectrality 
in the constant mass scenario has not yet been studied. It will be worthwhile to discuss 
physical acceptability of such an isospectrality in the position dependent mass background. 
In some articles (9 12 



21 
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61], solution of the PDM Schrodinger equation has 



been obtained without transforming the potential in to mass space . In this case the energy 
spectrum of the PDM Hamiltonians are not isospectral with the constant mass Hamiltoni- 
ans. Therefore it is reasonable to develop a method for solving PDM Hamiltonian without 
transforming the potential into the mass space. 

However, the fundamental question remains open: how the potential is affected when it 
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is expressed in the mass space? To answer this question, one has to obtain a solution for the 
Schrodinger equation without transforming the potential to the mass space. In this article, 
we will obtain a semi-analytical solution of the Schrodinger equation without transforming 
the potential to the mass space. This is another reason, to build a realistic model for solving 
PDM Hamiltonian. 

It is well known that the study of the same mathematical problems from different point 
of view lead to the progress of the science and includes a lot of mathematical tastes. A 
technique based on asymptotic expansion of Taylor series |64j has recently been suggested 
to obtain eigenvalues of Schrodinger equation which improves both analytical and numerical 
determination of the eigenvalues. Asymptotic Taylor Expansion Method (ATEM) is very 
efficient to obtain eigenvalues of the Schrodinger equation because of their simplicity and 
low round off error. This method has been easily applied to establish eigenvalues and wave 
function of the Schrodinger type equations. We would like to mention here that the ATEM 
is a field of tremendous scope and has an almost unlimited opportunity, for its applications 
in the solution of the constant and PDM Schrodinger equations. In this paper, we address 
ourselves to the solution of the PDM Schrodinger equation by using the ATEM. 

The paper is organized as follows. In the next section we review construction of ATEM by 
reformulating the well known Taylor series expansion of a function that satisfies second order 
homogeneous differential equation of the form: f"(x) = p (x)f'(x) + q (x)f(x). Section 3 is 
devoted to the application of the main result for solving the PDM Schrodinger equation for 
various forms of the Kinetic energy operator. As a practical example, we illustrate solution 
of the PDM Schrodin ger equation including harmonic oscillator potential and variable mass 
m{x) = mo (1 + 7a; 2 ) 11]. In this section, we present an approximate analytical expression 
for eigenfunction and numerical results for eigenvalues of the PDM Schrodinger equation. 
We also analyze the asymptotic behavior of the Hamiltonian. Some concluding remarks are 
given in section 4. 



FORMALISM OF ATEM 



In this section, we show the solution of the Schrodinger type equation for a quite ample 
class of potentials, by modifying Taylor series expansion by means of a finite sequence instead 
of an infinite sequence and its termination possessing the property of quantum mechanical 
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wave function. Let us consider Taylor series expansion [66] of a function f{x) about the 
point a: 

f(x) = f(a) + (x- a)f'(a) + \{x - aff{a) + \{x - a)*f®(a) + 

2 6 

CO / \ n 

n=0 

where f( n \a) is the n th derivative of the function at a. Taylor series specifies the value of 
a function at one point, x, in terms of the value of the function and its derivatives at a 
reference point a. Expansion of the function f(x) about the origin (a = 0), is known as 
Maclaurin's series and it is given by, 

f(x) = / ( ) + x/(0) + ix 2 r(0) + ^ 3 / (3) (0) + 



= E^/ (n) (°)- ( 2 ) 

n=0 

Here we develop a method to solve a second order linear differential equation of the form: 

f"(x)=p (x)f(x) + q (x)f(x). (3) 

It is obvious that the higher order derivatives of the f(x) can be obtained in terms of the 
f(x) and f'(x) by differentiating fl3]). Then, higher order derivatives of f(x) are given by 

f in+2 \x)= Pn (x)f(x) + q n (x)f(x) (4) 

where 

p n (x) = p (x)p n - 1 (x) +p' n _i(x) + q n -i(x) , and 

q n (x) = q (x)p n „x(x) + q' n _ 1 (x). (5) 

Of course, the last result shows there exist a formal relation between asymptotic iteration 
method (AIM) 7(1 and ATEM. We have observed that eigenfunction of the Schrodinger 
type equations can efficiently be determined by using ATEM. It is clear that the recurrence 
relations (jSJ) allow us algebraic exact or approximate analytical expression for the solution 
of ([3]) under some certain conditions. Let us substitute (jSJ) into the ([T|) to obtain 

(m n \ / m n \ 

1 + E^(0)^ + /'(°) 1 + E^(°)^[ • ( 6 ) 
n=2 ' / V n=2 ' / 



After all, we have obtained useful formalism of the Taylor expansion method. In the solution 
of the eigenvalue problems, truncation of the the asymptotic expansion to a finite number of 
terms is useful. If the series optimally truncated at the smallest term then the asymptotic 
expansion of series is known as superasymptotic (ti|, and it leads to the determination 
of eigenvalues with minimum error. Then boundary conditions can be applied as follows. 
When only odd or even power of x collected as coefficients of /(0) or /'(0) and vice verse, 
the series is truncated at n — m then an immediate practical consequence of these condition 
for g m _ 2 (0) = or p m _2(0) = 0. In this way, one of the parameter in the g m _2(0) and/or 
Pm-2(0) belongs to the spectrum of the Schrodinger equation. Therefore eigenfunction of 
the equation becomes a polynomial of degree m. Otherwise the spectrum of the system 
can be obtained as follows: In a quantum mechanical system eigenfunction of the system is 
discrete. Therefore in order to terminate the eigenfunction f(x) we can concisely write that 

9m(0)/(0)+p m (0)/ , (0) = 
g m _ 1 (0)/(0)+p m _ 1 (0)/ , (0) = (7) 

eliminating /(0) and /'(0) we obtain 

9m(0)p ro _i(0) - p m (0)g m _!(0) = (8) 

again one of the parameter in the equation related to the eigenvalues of the problem. 

In quantum mechanics bound state energy of the atom is quantized and eigenvalues 
are discrete and for each eigenvalues there exist one or more eigenfunctions. When we 
are dealing with the solution of the Schrodinger equation, we are mainly interested in the 
discrete eigenvalues of the problem. The first main result of this conclusion gives necessary 
and sufficient conditions for the termination of the Taylor series expansion of the wave 
function. 

The process presented here is iterative and number of iteration is given by m. The 
results are obtained as follows: in our Mathematica program, we use an iteration number, 
say m = 30, then we obtain another result for m = 40, so on, then we compare values of 
the parameter (eigenvalue) in each case till 10 digits. If values of the parameter reach its 
asymptotic value then we use these values and omit the others. For instance, if one can 
obtain values of the parameters for m = 40, first few of them will be reached its asymptotic 
values, say first 8 values. The following comment for the function is considerable: for such 
a solution it is suitable to take sum of first 8 term in the (jSJ). 
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It will be shown that ATEM gives accurate results for PDM Schrodinger equations. In 
the following sections, it is shown that this approach opens the way to the treatment of 
PDM Schrodinger equation including large class of potentials of practical interest. 



SOLUTION OF THE PDM SCHRODINGER EQUATION BY USING ATEM 

In the PDM Schrodinger equation the mass and momentum operator no longer commute, 
so there are several ways to define kinetic energy operator. The general expression for the 
Hamiltonian with the kinetic energy operator introduced by von Roos 
energy V(x), can be written as: 



67] and potential 



H = - (m n pm £ pm p + m p pm e pm v ) + V(x) (9) 

where r] + e + p = — lisa constraint and m = m(x) is position dependent mass. There are 
many debates for the choice of the parameters 77, e,and p, in our approach, we will obtain 



the solution of the PDM Schrodinger equation for the following Hamiltonians 



67J469J: 



Hi = \ ( P-P J + V(x); for e = -1, p = 0, 77 = 0, (10a) 
2 \ m J 

H 2 = U-p 2 + p 2 -)+V(x); fore = 0,p = 0,77=-l, (10b) 
4 \m m J 

H 3 = If-Lpa-Lj+y^); for £ = -i,p = 0,77 = -^ (10c) 



Hi = I ( p i p i + i p i p ) + nx) ' for £ = °< " = 4 " = 4 (10d) 



Here we take a new look at the solution of the the PDM Schrodinger equation by using the 
method of ATEM developed in the previous section. 

Before going further we share one of our significant observation during our calculations. 
If the mass distribution is not appropriate for a given potential, the eigenvalues do not reach 
their asymptotic values and resultant eigenfunction cannot be terminated when x — > ±00. 
In order to illustrate semi analytical solution of the eigenvalue equations 

H^(x) = Eij(x),(t = 1,2,3,4) (11) 

including harmonic oscillator potential, V(x) = ^m uj 2 x 2 , we use the mass distributions 
m(x) = m (1 + 7X 2 ), where 7 is arbitrary positive constant. By the way, we emphasize 
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that the wave function of harmonic oscillator potential is well defined in the region of ±00 

and satisfy that lim > 0. In this limit the mass distributions to be continuous. 

It is well known that asymptotic behavior of constant mass Schrodinger equation including 
harmonic oscillator potential is given by tp = e~~ f(x), for simplicity we set h = m — oj — 

l.Thus, this change of wave function guaranties lim y 0. After this transformation, 

we present an iteration algorithm to calculate both eigenvalues and eigenfunctions of the 
eigenvalue equation ffTTj) . Using this algorithm, we develop a Mathematica program, which 
demonstrates that it is easier to be implemented into a computer program, and produces a 
highly accurate solution with analytical expression efficiently 



Asymptotic Analysis 

The term asymptotic means the function approaching to a given value as the iteration 
number tends to infinity. By the aid of a Mathematica program we calculate eigenvalues 
and eigenfunction of Hi for 7 = 0.1 using number of iterations k = {20,30,40,50,60}. The 
function f(x) for n = 2 state is given in (j!2p and eigenvalues are presented in Table I. 



k = 20; f(x) = 1 - 1.857a: 2 - 1.619 x lO^V + 2.060 x 10"V 

+1.515 x HT 3 a: 8 - 1.261 x 10~V° - 6.495 x HTV 2 

k = 40; f(x) = 1 - 1.856x 2 - 1.622 x lO^V + 2.051 x 10~V 

+ 1.505 x KT 3 x 8 - 1.271 x 10 _4 x 10 - 6.662 x KTV 2 (12) 

k = 60; f(x) = 1 - 1.856x 2 - 1.622 x 10~ V + 2.051 x 10~V 

+ 1.504 x 10~V - 1.271 x 10"V° - 6.663 x 10~V 2 

Our calculation gives an accurate result for first 8 eigenvalues and eigenfunctions after 40 
iterations. Here we have used 60 iterations. Figure 1 shows the plot of normalized wave 
functions for first 6 state. 



Solution of the Hamiltonians H2 , H3 and H4 



In the previous section we have illustrated applicability of our method by solving Hamil- 
tonian H\. In this section we apply the same procedure to solve the Hamiltonians H2, H3 
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k 


n = 


n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


20 


0.46889047 


1.43341211 


2.35765542 


3.28397486 


4.21360362 


4.35399596 


30 


0.46889665 


1.43348058 


2.35642259 


3.24660834 


4.12086916 


4.98321327 


40 


0.46889650 


1.43348582 


2.35655507 


3.24585555 


4.10543833 


4.95755341 


50 


0.46889651 


1.43348553 


2.35654885 


3.24599291 


4.10703835 


4.94114551 


60 


0.46889651 


1.43348555 


2.35654908 


3.24598255 


4.10694346 


4.94337909 



TABLE I: Eigenvalues of the PDM Hi for different iteration numbers k and 7 = 0.1. 

and H4. Again we have used 60 iterations for each Hamiltonians and checked stability of the 
eigenvalues. Here we calculated eigenvalues for 30 iterations and they are listed in Table II. 
We have also checked that for the given eigenvalues, the wave functions are normalizable 
and it tends to zero when x — > 00. 





E H 2 


Eh 3 


Eh 4 





0.50773226 


0.48833347 


0.50949336 


1 


1.45551369 


1.44451856 


1.45972923 


2 


2.36941282 


2.36286881 


2.37461896 


3 


3.25544187 


3.25137213 


3.26106459 


4 


4.13235379 


4.12882619 


4.13805287 


5 


4.95997506 


4.96305356 


4.96478901 



TABLE II: The eigenvalues Hamiltonians H2, H3 and H^, for 7 = 0.1. The result is obtained after 
30 iterations. 

The results given in Table II shows that eigenvalues and eigenfunctions are also depends 
on the choices of the parameters, e, p, and 77 of Hamiltonian OU). 

REMARKS AND DISCUSSIONS 

In this paper, we have studied the solution of the PDM Schrodinger equation without 
mapping the potential in to the mass space. We have solved PDM Schrodinger equation 
for four different kinetic energy operators including harmonic oscillator potential with the 
variable mass function of the form m(x) = mo (I + *~yx 2 ). It is shown that energy levels 
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FIG. 1: Plot of the normalized wave function of the PDM Hamiltonian (|10a|) for n = 0, 1, 2, 3, 4, 5. 



of the PDM Schrodinger equation depends on the mass distributions. It is important to 
remark that the results presented here, shows that eigenvalues also depends on the ordering 
parameters of the PDM Schrodinger equation 

We have presented an approximate method based on asymptotic Taylor Series Expansion 
of a function. Fortunately, this method is useful for obtaining both eigenvalues and eigen- 
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functions of the Schrodinger type equations. Therefore, the results have been obtained here, 
allowing further comparisons between the models. 

As a further work the method presented here can be used to built more realistic models 
for the PDM physical systems. Before ending this work a remark is in order. When the 
potential mapped to the mass space, the both constant and PDM Hamiltonian has the same 
eigenvalues. It will be worthwhile to discuss physical acceptability of such an isospectrality 
in the position dependent mass background. Therefore we have to develop methods for 
solving PDM Schrodinger equation without connecting mass to potential or vice versa. 
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